Hadronic contribution to the muon g-2: a theoretical determination 
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The leading order hadronic contribution to the muon g-2, a^ AD , is determined entirely from theory 
using an approach based on Cauchy's theorem in the complex squared energy s-plane. This is 
possible after fitting the integration kernel in a^ AD with a simpler function of s. The integral 
determining a^ AD in the light-quark region is then split into a low energy and a high energy part, 
the latter given by perturbative QCD (PQCD). The low energy integral involving the fit function 
to the integration kernel is determined by derivatives of the vector correlator at the origin, plus a 
contour integral around a circle calculable in PQCD. These derivatives are calculated using hadronic 
models in the light-quark sector. A similar procedure is used in the heavy-quark sector, except that 
now everything is calculable in PQCD, thus becoming the first entirely theoretical calculation of 
this contribution. Using the dual resonance model realization of Large N c QCD to compute the 
derivatives of the correlator leads to agreement with the experimental value of a M . Accuracy, though, 
is currently limited by the model dependent calculation of derivatives of the vector correlator at the 
origin. Future improvements should come from more accurate chiral perturbation theory and/or 
lattice QCD information on these derivatives, allowing for this method to be used to determine 
a^ AD accurately entirely from theory, independently of any hadronic model. 
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The value of the muon g-2 is well known as a test of the 
standard model (SM) of particle physics. The SM 
result for the anomalous magnetic moment of the muon 
is conveniently separated into the contributions due to 
QED, the hadronic sector, and the electroweak sector. 
A sizable theoretical uncertainty arises from the (leading 
order) hadronic vacuum polarization term, the second 
largest contribution after that of QED. A substantial ef- 
fort has been made to determine this contribution from 
experimental data on e + e~ — > hadrons and r — » hadrons 
@-[||. Currently, there is some yet unresolved discrep- 
ancy between both results. Writing the muon anomaly 
in the SM as 



a S , M = a* ED 



HAD 



(1) 



the leading contribution is that from QED, followed by 
the hadronic and the electroweak terms. In this paper 
we concentrate on the leading order hadronic contribu- 
tion and discuss a new approach to its calculation entirely 
from theory. The method relies on Cauchy's theorem in 
the complex squared energy s-plane, after fitting the in- 
tegration kernel entering a^ AD with a simple function 
of s . In the region of the light-quark sector the method 
requires knowledge of some of the derivatives of the (elec- 
tromagnetic) vector correlator at zero momentum, as well 
as its perturbative QCD (PQCD) behavior. Currently, 
these derivatives will be obtained here from hadronic 
models, thus being affected by systematic uncertainties. 
Hence, at this stage the method cannot rival in accu- 
racy with the standard approach of using experimental 
data on the vector correlator at low/intermediate ener- 
gies. However, future precision determinations of these 
derivatives from chiral perturbation theory (CHPT) or 
lattice QCD would render this calculation of a^ AD inde- 
pendent of experimental data on e + e~ —¥ hadrons and 



t — > hadrons. It must be emphasized that as a conse- 
quence of Cauchy's theorem this method only requires 
knowledge of the derivatives of the vector correlator at 
the origin, rather than its full expression over an ex- 
tended energy region. In addition, the method allows 
for a straightforward incorporation of the charm- and 
bottom-quark contributions to calculable exclu- 

sively from PQCD, i.e. without the need for data on the 
vector correlator. This leads to the first entirely theoret- 
ical calculation of this contribution. 
We begin with the standard expression of the (lowest or- 
der) muon anomaly [l[ 



HAD 



3tt 2 



ds 

™ K{s) R{s) 



(2) 



where olem is the electromagnetic coupling and the stan- 
dard R -ratio is R(s) = 3 J2f Q} [8 n ImII(s)] , where Q / 
are the quark charges and n(s) is the vector current cor- 
relator normalized to 8 7rlmn(s) = 1 + a s /n + • • • . The 
integration kernel K(s) in Eq.© is given by [H 



K(s) 



dx ■ 



x 2 {l-x) 



-(l-x) 



(3) 



where m M is the muon mass. A popular approach to 
compute a^ AD has been to split the integral in Eq.((2j 
into a low energy region from threshold up to s — sq — 
(1.8 GeV) 2 , followed by a high energy region from s = sq 
to infinity. The integral in the former region was cal- 
culated using data on of e + e~ — > hadrons or (isospin 
rotated) data on r — > hadrons. The integration in the 
high energy region was performed assuming PQCD. In 
this paper we discuss a new approach based entirely on 
theoretical input. For later convenience we split the con- 
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tributions to the leading order a^ AD into three pieces, 
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HAD I 



\uds 



HAD \ 



HAD I 
a u \b 



(4) 



where the first term on the right hand side corresponds 
to the contribution of the three light quarks, and the 
second and third term refer to the charm- and bottom- 
quark contributions. In the light-quark sector the first 
step is to fit the integration kernel K (s) in an interval 
sth < s < sq with a function K\(s) 



Kx{s) = a s + 



Ea r , 



(5) 



with coefficients determined by minimizing the chi- 
squared. The upper limit so is below the charm thresh- 
old. Next, the integration range in Eq.© is split into 
a low energy (s < sq) and a high energy (s > sq) re- 
gion where PQCD would be valid. In the former region 
Cauchy's theorem is used to obtain 



f s ° ds 1 

/ — -ImlWa) =Res 



Ki(s) 



1 f ds 

— <b — Ki(s) Tluds(s) 
2m J\ S \ =SQ s 



s=0 

(6) 



where the integral on the right hand side, around the cir- 
cle of radius so ~ (1.8 GeV) 2 , is computed using PQCD 
in the light quark sector. This is known up to five-loop 
level [j| . The contour integration can be performed using 
fixed order perturbation theory (FOPT) or, alternatively, 
contour improved perturbation theory (CIPT). There is 
no clear a priori criterion to decide which is best in a given 
application. However, in the present case the difference 
between the two methods turns out to be negligible, as 
discussed later. The residues are given in terms of deriva- 
tives of the correlator at zero momentum, for which one 
can use hadronic models, CHPT or lattice QCD. Hence 
Eq.© becomes 



n HAD\ 
u p, \uds 



t—u,d,s 

— I *L 

27TJ /i.i_.„ S 



n„<j s (s) 



s=0 



Ki(s) n uds (s)\pQCD 



J \ s\— s 

ds 1 
— K(s) -Im U uds (s) | pqcd\, (7) 



where the last integral above involves the exact inte- 
gration kernel K(s) and PQCD is used for the spectral 
function. It is important to stress that this contribu- 
tion to the anomaly only requires knowledge of a few 
derivatives of the vector correlator at the origin (to 
compute the residue). It does not require knowledge of 
the correlator itself in the extended energy region from 
threshold up to s — (1.8 GeV) 2 . The choice of this 
particular value for the onset of PQCD will allow for a 
fair comparison with determinations based entirely on 
data 0. It is also supported by experimental results 



from the BES Collaboration 6] which show the onset of 
PQCD at so ~ 4.0 GeV 2 . The stability of results against 
changes in this threshold value will be analyzed later. 



In order to incorporate charm-quark information we add 
an extra contribution determined as follows. A new fit 
to the integration kernel K{s) is performed in a region 
si < s < S2, where s\ ~ Mj,^, and s 2 ~ (5.0 GeV) 2 . 
Using this kernel and Cauchy's theorem the charm con- 
tribution is given by an expression similar to Eq.©, ex- 
cept that so is replaced by s 2 and Ki(s) by the new fit 
function ^(s). The residues can now be computed di- 
rectly from PQCD using the low energy expansion of the 
heavy-quark correlator, known up to four-loop order. No 
hadronic model nor data is needed here. A similar pro- 
cedure can be followed to incorporate the contribution of 
the bottom quark. 

We proceed to fit the integration kernel, Ki(s) in the re- 
gion sth < s < so- If one were to choose a polynomial fit 
of the form Ki(s) = J2i=i c i s% then the residues in Eq.© 
would all vanish and the anomaly would be determined 
entirely from QCD Q. There are two drawbacks to such 
a fit. First, even taking many terms in the series expan- 
sion of Ki(s) the fit is not accurate enough. Second, the 
higher powers of s bring in higher dimensional conden- 
sates in the operator product expansion, thus reducing 
further the accuracy of this approach. An inspection of 
the s-behavior of the kernel K(s) suggests that a series 
expansion involving inverse powers of s should be a bet- 
ter option. In fact, this turns out to be the case, e.g. the 
fit function, Eq.©, up to s~ 3 becomes 



ifi(s) = 2.257 x 10~ 5 s + 3.482 x lO^s" 1 

- 1.467 x 10~ 4 s~ 2 + 4.722 x lO^s" 3 , (8) 

where s is expressed in GeV 2 , and the numerical coef- 
ficients have the appropriate units to render K\(s) di- 
mensionless. Figure 1 shows the exact kernel K[s) in 
Eq-© (solid curve) together with the fit -f^i(s) as in 
Eq.© (solid dots). The relative difference between the 
two curves lies in the range — 1% in the low energy re- 
gion, where it contributes the most. A further estimate 
of the accuracy of the fit function, Eq.©, can be ob- 
tained by using all available experimental data on R(s) 
in Eq.© together with (a) the exact kernel Eq.©, and 
(b) the fit kernel Eq.©. We find af[ AD \ uds = 641.69 



for procedure (a) and a 



HAD I 



\uds 



641.16 for procedure 



(b), i.e. a difference of 0.08%. Using additional inverse 
powers of s terms in the fit, while improving it slightly, 
it does not lead to any appreciable difference in the final 
result for a^ AD . For instance, the difference in a^ AD 
from adding two additional inverse powers in Eq.© is 
less than 0.16%. 

We discuss now the incorporation of the heavy quark 
sector contribution to the anomaly, starting with the 
charm-quark piece a^ AD \ c . The fit to the integration 
kernel, which we name ^(s), is performed in the region 
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differs from the exact kernel by less than 0.0005 % in the 
range My < s < (12 GeV) 2 . The residue is now 



Res 



H-b(s)\pQCD 



K 3 (s) 



= 6.3 x 10" 



(14) 



1.5 2.0 
s (GeV 2 ) 



where the error is negligible. Next, in order to calculate 
the contour integral around the circle of radius s 2 we 
make use of PQCD, i.e. 

np Q cD( s ) = E(^) n(n) W, (is) 

n=0 ^ ' 



where 



FIG. 1: The exact kernel K(s), Eq.^jfy (solid line) together with 
the fit in the light-quark region, Eq. (JSJ, (solid circles). 
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function 



< s < s 2 



K 2 (s) 



(5.0 GeV) 2 . The very simple 



s 



«2 
e 2 



(9) 



where a t = 0.003712 GeV and a 2 = -0.0005122 GeV 4 , 
provides an excellent fit with K 2 (s) differing from the ex- 
act kernel K(s) by less than 0.02%. In this case a^ AD \ c is 



given by a similar expression as a 



HAD I 



\uds i 



Eq.©, with 



The complete analytical result in PQCD up to 
0(a 2 , (m 2 /s) 6 ) is given in [13], with new results up to 
order 0(a 2 (m 2 / s) 30 ) obtained recently [l5j]. There are 

also exact results for li^ and I Ii ' from [lj| , while Tl 2 3 ^ 
is known up to a constant term [17j ] . This constant term 
does not contribute to the contour integral due to the 
s-dependence of K 2 {s). Finally, at five-loop level the full 
logarithmic terms in Tl^ and nj 4 ' are known from [l8[ 



obvious replacements. An important difference is that and respectively. The contour integrals in FOPT 



,HAD\ 



only involves the correlator and its imagi- 
nary part calculable in PQCD, i.e. no hadronic model is 
needed for the residue. In fact, the Taylor expansion of 
the correlator around the origin is given by 



1 f ds 

f ~ Kn ^ U i( s )\pQCD 



135.3(6) x 10~ 7 
20.3(1) x 10~ 7 
3.6(2) x 10~ 7 



(17) 



Kc{s)\pQCD 



32ir- 



Q 



n>0 



C n z n 



(10) 



where z = s/(4m 2 ). Here m c = m c (/i) is the charm- 
quark mass in the MS-scheme at a renormalization scale 
fi. The coefficients C n up to n = 30 are known at 
three-loop level Q-B. At four- loop level C and C\ 
were determined in jgj-^3|, C 2 in [9 1 ] and C 3 in fTTj j . 
Due to the s-dependence of K 2 (s) no coefficients C4 and 
higher contribute to Res[n c (s)p(s), s = 0]. Using as 

input fi = 3 GeV, ai 4) (3GeV) = 0.2145(22) [H and 
m c (3 GeV) = 0.986(10) GeV we find 



for n — 1,2,3 and q = uds,c,b, respectively. For n= 1 
the result in CIPT is 135.6(6) x 10" 7 , i.e. a 0.2% dif- 
ference with FOPT. Also for n = 1, changing the PQCD 
threshold in the interval sq = (1.8 — 2.0) 2 GeV 2 leads to 
a change of only 0.15% in the final value of a^ AD . The 
BES Collaboration data [6J agrees well with PQCD in 
this region and beyond. The results for the third inte- 
gral in Eq.([7]), and their equivalent for the charm- and 
bottom-quark sectors are 



f°° ds , s l 

/ — K(s)-lmn q {s)\ PQC D 

J Si s 7T 



n c (s) 



0.03604 + 0.001833 s + 0.00012335 s 2 
0.000012472 s 3 + 0{s 4 ) , (11) 



151.8(1) x 10- 7 
20.0(4) x 10" 7 (18) 
3.4(2) x 10~ 7 



where s is expressed in GeV , and the numerical coeffi- 
cients have the appropriate units to render II c (s) dimen- 
sionless. The residue in the charm-quark sector is 



with j = 0, 2, 4 for q = uds, c, b, respectively. Substitut- 
ing the results from Eas.([I7jl and (IT5|) into Eq.JT]) and the 
corresponding expressions for the charm- and bottom- 
quark contributions, the leading order a 



HAD 



Res 



n c (s)|pQ C D 



K 2 (s) 



76.1(5) x 10~ 7 , (12) 



HAD 



s=0 



where the error is due to the uncertainty in a s and to 
the truncation of PQCD. For the bottom quark case the 
fit function 

K 3 (s) = 0.003719 GeV 2 s' 1 - 0.0007637 GeV 4 s~ 2 (13) 
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= y a iM Res n ud;5 (s) 
+ 19.4(2) x 10~ 10 . 



KM 



(19) 



The contributions to a^ AD from the charm- and bottom- 
quark sectors obtained from PQCD are 

a" AD \ c = 14.4(1) x 10- 10 , (20) 
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and 



a" AD \ b = 0.29(1) x 10" 



(21) 



Finally we discuss the calculation of the hrst term on the 
right hand side of Eq. ([6]) . Given the parametrization in 
Eq.© this term can be conveniently written as 



Res 



>o ^— ' n! ds r ' 

n—l 



where the 



are the coefficients of the s 



and 



s~ 3 terms in Eq.©, respectively. To demonstrate the 
effectiveness of the method we consider three hadronic 
models for the vector correlator, (single p) Vector Meson 
Dominance (VMD), the Kroll-Lee-Zumino (KLZ) quan- 
tum field theory model [20| -[21], and the dual reso- 
nance model realization of QCD in the large N c limit 
(Dual-QCD co )[22]-[2l|. We reiterate that the use of 
hadronic models to compute the derivatives of the vec- 
tor correlator at the origin is only provisional. In future, 
these derivatives will be provided with increased accu- 
racy by CHPT and/or lattice QCD. The error on VMD 
can be estimated to be of order O(10 — 20%) judging 
from its predictions of the pion radius and form factor. 
The KLZ model is a renormalizable theory of pions and 
a neutral p-meson which provides the necessary quantum 
held theory platform for VMD, and leads to loop correc- 
tions to VMD. The loop corrections in the KLZ model 
bring the pion radius and form factor into better agree- 
ment with experiment. In VMD-type models the vec- 
tor correlator is related to the pion form factor through 
IWs) = F„{s)/tf, where f p = 4.96 ± 0.02 is the 



7 — p coupling and F w (s) is the pion form factor. 
p-VMD expression for the correlator is 



Huds(s)\ 



VMD 



M 2 

±v± p 



1 



n {mi 



s) 



The 



(23) 



which involves the underlying standard VMD universal- 
ity relation qp^ j ' j " p = 1. With g p7r7r = 5.92 ± 0.01 from 
experiment [l2|, this relation is off by roughly 20%. The 
result for the residue in Eq.(fT9| is 



Res [n uds (s)K 1 (s)/s 
leading to 



VMD 



2.20(2) x 10" 



HAD\ 



\VMD 



644(6) x 10 



-10 



(24) 



(25) 



For the correlator in the KLZ model we use the result 
from (25| (see also [2JJ) and obtain 



Res [Iluds&K^/s}^ 2 = 2.22(2) x 10 



and 



HAD I 
a p \KLZ 



650(6) x 10 



-10 



(26) 



(27) 



The errors for VMD and KLZ are only those due to the 
uncertainty in f p and do not include possible (system- 
atic) model errors. The latter can be gauged from the 
deviation from universality gpmr/fp = L °ff by some 
20%, as well as from the pion charge radius in VMD 
< r 2 >= 0.394 fm 2 , to be compared with the experimen- 
tal value [HI < r 2 >= 0.439 ± 0.008 fm 2 . 
While QCD in the limit of an infinite number of col- 
ors leads to a hadronic spectrum consisting of an infi- 
nite number of zero-width resonances, it does not specify 
the mass spectrum nor the couplings. Dual-QCD^ [221 ]- 
[24| provides this information leading to hadronic form 
factors in excellent overall agreement with data in the 
space-like region. The vector correlator in this frame- 
work is given by 



n„ds(s)|Qcz5 oc - -pi 
J P 



i i r(/3-i/2) 



X 



tt r(/3-i) 

B(/3-l,l/2- s /2M 2 ) 



(28) 



where (3 is a free parameter and B(x,y) is the Euler beta- 
function. From the power series expansion of B(x,y) it 
is easy to see that Eq. (|28l) represents an infinite num- 
ber of (zero-width) resonances corresponding to the p- 
meson and its radial excitations. The latter account for 
the deviation from the VMD result g p7 m/ f P = 1 lead- 
ing to [22] <? p7r7r / fp = 1.2 in agreement with experiment. 
For p = 2 Eq.® reduces to single p-VMD. The value 
f3 = 2.30(3) results in an excellent fit to all data on 
the pion form factor F 7r (s) in the space-like region up to 
s = — lOGcV 2 with a chi-squared per degree of freedom 
Xf ^ 1-5 In contrast the VMD fit yields xf ^ H- 
In addition, the Dual-QCDoo model gives a value of the 
pion charge radius < r 2 >= 0.436 ± 0.004 fm 2 H3 to be 
compared with the most recent experimental value [26j 
< r 2 >= 0.439 ± 0.008 fm 2 . Since the first derivative 
of the vector correlator dominates in Eq. (|22|) . it is very 
important for a hadronic model to reproduce the pion 
radius. The result for the residue in Eq. (fT9| is 



Res 



QCD a 



= 2.47(3) x 10" 



and using Eq. ([T9l the hadronic af- CDoc is 



HAD 



\ QCDoo = 722 (9) x 10 



-10 



(29) 



(30) 



where the error is mostly due to that in j3. The 
result Eq. ([3H|) can be compared with the value 
af[ AD = 692.3(4.2) x 10~ 10 from 0-0 using e+e" 
data, or a" AD = 701.5(4.7) x lO" 10 using r-data. How- 
ever, a more recent reanalysis based on r-data [27T ] finds 
a^ AD = 690.96(4.65) x 10~ 10 . In the QCD^ framework 
the 1/N C corrections arise in the time- like region from 
finite width resonance effects. These corrections to the 
form factor are of order C(r 2 /M 2 ) near the origin. 
While small, they might have an impact on the residues, 
Eq.flgS]). Shifting the poles in Eq.® to the second 
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Riemann sheet in the complex s-plane, while preserving 
the normalization at the origin, and the vanishing of the 
imaginary part of F 7r (s) at threshold [11], [11] leads to a 
reasonable finite width model. The first derivative of the 
form factor at the origin, which is the main contribution 
to Eq. ([29|) . receives no width correction. The second 
derivative is reduced with respect to the zero-width 
result by less than 2%, and the third derivative by some 
3%. This translates into an increase in the value given 
in Eq.® of 0.1%. 



Adding to Eq.flUJl the QED contribution [H a^ ED = 
11658471.809±0.015, the electro-weak [U a% w = 15.4± 
0.2, the higher order hadronic a^ AD \ HO = -9.79 ± 
0.09, and the light-by-light contribution [31] a^ bL \ HO = 
11.6 ±4.0, all in units of 10~ 10 , we find it intriguing that 
the Dual-QCDoo prediction, EqJ5Dl leads to 



leading to 



HAD 



815 x 10 



-10 



(35) 



This is a great deal larger than the value expected from 
experiment. The reason for this is that the constant Cg 3 
dominates this result. Furthermore, in [32| it is argued 
that the estimates of the 0(p 6 ) constants could be larger 
than the physical constants. Therefore it makes sense 
that this result for a^ AD represents an overestimate, 
rather than an underestimate. One can reverse this 
argument and give the first model independent determi- 
nation of C% 3 . Making use of af[ AD = 692.3(4.2) from 
[2], we find that C£ 3 = -13.9(2) x lO^GeV" 2 . 



a^QCD^ = 11659210.6 ±9.8 x 10" 1 
to be compared with the experimental value 



a EXP = 11659208.9 ±6.3 x 10" 



in 



Equation (f3ij) does suggest that it might be possible to 
understand the muon anomaly entirely within the SM. 
Our approach to determine a^ AD appears to be opti- 
mally designed for use in CHPT, as the main input is 
the power series of the correlator around the orig in. The 
Q(p 6 ) vector correlator was determined in [32| and in 
[33j. The derivative at zero momentum, in terms of the 
usual chiral constants, is 



as 



IxpT 



(0) 



0.0105557- 
1.0346LS . 



ACL 



0.77725L^ 

(33) 



Two of these constants have been calculated on the 
lattice, LI = 3.08(23) (51) x 10" 3 ja| and L r w = 
-5.2(2) (if) x 10~ 3 [35]. The constant C£ 3 has not been 
determined on the lattice yet. We have to rely on a very 
rough estimate of this constant from [32J using VMD, 
C£ 3 » -17 X 10~ 3 GcV~ 2 . With these values, we obtain 
■j^U. u dsifi) | xpX ~ 0.0857. Given the very large uncertainty 
in Cg 3 , and the fact that the contribution of the second 
and the third derivative is very small, we find 



(31) In summary, we have discussed a new approach to the 
determination of the leading a^ AD entirely from theory, 
i.e. without the use of experimental data on the vector 
correlator in an extended energy region. This can be 

(32) achieved by fitting the integration kernel, Eq.Q, in the 
light-quark sector with the simple function Eq.((8]), and 
subsequently invoking Cauchy's theorem in the complex 
s-plane. This leads to the result Eq.0 which only 
requires knowledge of a few derivatives of the vector 
correlator at the origin. This must be contrasted with 
the standard approach which requires the complete 
correlator in the wide energy region from threshold up 
to so — (1.8GeV) 2 . Such a detailed information can 
only be reliably and accurately obtained from data. 
Currently, these derivatives can be estimated using 
hadronic models, examples of which have been presented 
here. In future, though, more accurate determinations of 
the derivatives from CHPT and/or lattice QCD should 
become available thus allowing for a model independent 
calculation of this contribution. In the heavy-quark 
sector this problem does not arise, as it is possible to 
calculate the anomaly entirely from PQCD, with the 
results given in Eas.(P0" ]) -([2"T |) . 



Res 



n„ds(s) 



s=0 



2.8, 
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